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Abstract 

If current ideas about unified field theories are correct, macroscopic cosmic 
defects may well exist. The observation of such an entity would have enormous 
significance for our understanding of fundamental physics. This paper points 
out a novel observable signature of cosmic texture and global monopoles, 
namely strong hot spots in the cosmic microwave anisotropy pattern on sub- 
degree scales. This signal should be readily detectable by the next generation 
of anisotropy mapping experiments. The signature arises from overdensi- 
ties in the photon-baryon fluid generated by the gravitational attraction of 
the defects. The angular power spectrum of the anisotropy fluctuations on 
subdegree scales is also calculated, for cosmic string, global monopoles, and 
texture. 

The prospect of accurate maps of the temperature fluctuations on the cosmic microwave 
sky gives us hope of resolving one of nature's most profound questions, the origin of large 
scale structure in the universe. The recently approved COBRAS-SAMBA and MAP satel- 
lites, the Very Small Array interferometer, and other experiments, promise to provide cos- 



mological data of unparallelled quality and quantity. 

Cosmic defects have been of fascination to theorists ever since it was realised that their 
presence was generic in unified field theories of high energy physics jl] . Indeed, if current 
ideas about the unification of forces are correct, cosmic defects in one form or another are 
almost unavoidable. Two such defects occur in standard model physics - namely electroweak 
textures and skyrmions, formed as a result of the breakdown of the electroweak and strong 
chiral symmetries respectively. Neither of these is relevant to the problem of structure for- 
mation, but analogous defects formed through the breakdown of symmetries at much higher, 
grand unification energy scales, do give rise to interesting theories of structure formation in 
the universe 0. Cosmic strings are the best known of these, but global monopoles and tex- 
tures are also promising candidates ||. A period of inflation can, but does not necessarily, 
dilute the defects to unobservably low densities. For example, recent work has has shown 
that they may be efficiently produced during the reheating process following inflation [|J. 
The imminent maps of the microwave anisotropy offer a unique observational window on 
cosmic defects, and the possibility of discovering whether any of these relics of unification 
have survived within our observable universe. 

As mentioned, macroscopic cosmic defects provide a possible mechanism for structure 
formation which is a valuable alternative to the more popular scheme based on inflation. 
The two theories could hardly be more different - where the inflationary scheme is inherently 
quantum mechanical, the cosmic defect mechanism is purely classical. The formation of 
fluctuations during inflation involves linear physics, whereas the defect theories are highly 
nonlinear. Related to this, the inflationary fluctuations take the form of Gaussian random 
noise, whereas the defect induced fluctuations are highly non-Gaussian. 

Recent work has emphasised further fundamental differences. First, the defect theories 
are highly constrained by causality, so that the fluctuations are completely uncorrelated 
beyond the horizon scale at all times |J, ||, ||. And second, the nonlinearity of the 
theories means that the evolution of each Fourier mode k in the perturbations is different, 
and cannot be predicted on the basis of the initial conditions of that mode alone, because 
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all the modes are coupled. This property is referred to as 'incoherence' - modes of the same 
wavenumber k do not follow identical time evolution J7[. This has the effect of smearing 
the 'Doppler' peaks in the angular power spectrum of the anisotropies, caused acoustic 
oscillations in the photon-baryon fluid just prior to recombination. 

Advances in computer power have recently made possible realistic calculations of mi- 
crowave anisotropy maps in the defect theories. Following the COBE satellite results, several 
groups have reported calculations for the large angle cosmic microwave anisotropy patterns 
T0|1 , |TT|| , |fl^1 , Unfortunately, the information available in a map with 7 degree res- 



olution, such as the COBE 4 year maps, appears too limited to meaningfully discriminate 
between the inflationary and defect theories - all these theories are compatible with the 
COBE measurements. 

The first realistic computations of degree-scale CMB anisotropy from cosmic defects 



were reported some time ago ||14|| . These were performed assuming a fully ionised universe, 
mainly for technical simplicity. The present paper shall consider instead the case of standard 
recombination, in which much more small scale structure is visible, on scales down to 10 arc 
minutes. Below this scale, even with standard recombination, photon diffusion on the last 
scattering surface begins to blur the anisotropy pattern. 

At early times, redshifts Z >> 1100, photons are tightly coupled to ionised matter, 
forming the photon-baryon-electron fluid. As the universe cools, the electrons and baryons 
undergo 'recombination' forming neutral atoms, and the plasma quickly becomes trans- 
parent. In the approximation that recombination happens instantaneously, the microwave 
anisotropy seen in a direction n on the sky is given by 

^"( n ) = - n • v 7 (i) - ^J f drh^M (1) 

where r is conformal time, 5 7 is the photon fluid density contrast, and v 7 the velocity of the 
photon-baryon fluid on the spherical shell representing the 'surface of last scattering'. The 
term |<5 7 is the intrinsic temperature contrast, and the second term is the Doppler effect. 
The last term involves the time derivative of the metric perturbation h^, and represents the 



change in the proper path length due to a time-dependent gravitational field along the line 
of sight. This is the Sachs- Wolfe effect §L§. 

Equation ([]]) is written in the synchronous gauge, sometimes criticised in the literature 
because it possesses residual gauge freedom. However, this gauge is well suited to the cos- 
mic defect theories discussed here. In these theories, the universe is assumed to be perfectly 
smooth before the symmetry breaking phase transition which produced the defects. The 
perturbation variables can therefore all be set to zero. This choice, called initially unper- 
turbed synchronous gauge has no residual gauge freedom. The nice property of this gauge 
is that the Einstein equations are then completely causal, with the value of all perturbation 
variables at a given spacetime point being completely determined by initial conditions within 
the point's backward light cone. Gauges like the Newtonian gauge suffer from the problem 
that perturbations due to anisotropic stresses, which are generally present in these theories, 
propagate acausally |16[| . 

In initially unperturbed synchronous gauge, all the variables in equation ([I]) have an 
unambiguous physical meaning. The contrast <5 7 in the photon density is that observed by 
a local freely falling observer (for example a cold dark matter particle). Likewise v 7 is the 
velocity of the photon-baryon fluid with respect to cold dark matter. There is an additional 
Doppler term which may be attributed to the peculiar motion of the dark matter, but that 
arises from an integration by parts of the last term [jl4}| . 

The split represented in equation (JJ) is useful in a causal theory, because the first two 
terms are completely determined by the source for perturbations within the backward sound 
cone, with comoving radius J Trec c s (r)<ir, with r rec the conformal time at recombination and 
c s < l/\/3 the speed of sound in the photon-baryon fluid. The last term depends on 
the perturbation source over a much larger region of spacetime, filling the backward light 
cone of each point on the photon's trajectory. Apart from the earliest contributions to the 
integral, in the vicinity of the last scattering surface, there is little correlation between the 
first two terms in (1) and the Sachs Wolfe integral. Furthermore, the latter produces an 
approximately scale invariant angular power spectrum, monotonically falling, and with no 
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'Doppler' peaks (see ref. |14|). It is therefore unlikely to be important in influencing the 
position of the peaks in the power spectrum produced from the first two terms, and is a 
subdominant effect in the height of the peaks. I have checked this explicitly in the texture 
model of reference ||, where the approximation of ignoring the Sachs- Wolfe integral in (1) 
yields a primary peak in the angular power spectrum higher than the full result, but only by 
20 per cent. Note that this approximation is certainly not good for a theory like inflation, 
where there are large correlations on superhorizon scales. 

In this paper I report on calculations of the first two terms alone, which dominate the 
anisotropy on scales from a degree down to several arc minutes. Including the Sachs- Wolfe 
integral is certainly a goal for future work, but the memory and CPU requirements would 
have severely compromised the resolution available. The calculational technique involves an 
improved nonlinear sigma model code in which local stress energy conservation holds to an 
accuracy of a few per cent, and is quadratically convergent in the timestep. This is coupled 
to a three dimensional leapfrog code for the dark matter (taken to be 'cold' dark matter), 
neutrinos (treated as a fluid with viscosity), and the photon-baryon-electron fluid, with a 
time-varying speed of sound. For strings, the flat spacetime code of ref. pi is employed, 



with the modification that an extra fluid (with c 2 s = 1/3) is introduced to account for the 
energy and momentum lost by the strings due to a) the expansion of the universe, and b) 
the emission of small loops and gravitational waves. This ensures exact local conservation 
of the source stress energy tensor. Whilst the string network used is reasonably realistic, 
having for example similar scaling density to that reported in expanding universe codes, the 
string results have considerably larger systematic uncertainties than those for texture and 
monopoles. 

Simulations were performed for 256 3 boxes, with the horizon at last scattering represent- 
ing 45 grid spacings. Tests of the source stress energy tensor indicate that the fluctuating 
part achieves scaling form after as little as 5 grid spacings, which is still significantly be- 
fore matter-radiation equality. At the instant of recombination, slices across each box were 
taken at intervals of 32 grid spacings (no correlation was apparent between neighbouring 



slices), with the Doppler and intrinsic terms stored separately. This produced 8 maps per 
simulation. The maps are available via ftp from the author. 

The four colour Figures show representative maps for strings, global monopoles, textures 
and standard cold dark matter. The maps are 10° square, and have been smoothed with a 
Gaussian of FWHM 12" (5 grid spacings), to simulate the effect of Silk damping. The maps 
show the temperature contrast (5T/T) in units of one standard deviation, arising from the 
first two terms of equation (1). In each case the cosmo logical parameters used were Q = 1, 
Hb = 0.05, h = 0.5, A = 0. As mentioned above, because the Sachs Wolfe integral has been 
omitted, these maps have slightly less power than they should on scales above one degree or 
so. But the most striking features in the maps occur on smaller scales, and these should be 
accurately represented. 

There are two ways in which the theories studied are manifestly different. First, the 
scale of the structure present in the maps increases as one goes from strings to monopoles, 
textures and standard inflation. This is a reflection of the differences in angular power 
spectra, and in particular the shift in the primary peak in the power spectrum (Figure 5). 
Second, the texture and monopole maps show clear 'hot spots' in excess of 5 or 6 a, which 
would be highly improbable in a Gaussian theory. The fact that hot spots, and no cold 
spots are seen, is a simple result of the fact that the defects themselves represent localised 
concentrations of energy, which attract the photon-baryon fluid, and cause a local rise in the 
local density and temperature. Localised defects, like unwinding textures or annihilating 
monopole-antimonopole pairs, do so more effectively than extended objects like strings. 
After a hot spot is produced, I have established that it is the presence of the somewhat 
diffused surviving defect stress energy in the region that prevents the photon density from 
oscillating negative and producing a cold spot. 

Let us now discuss quantitative measures of these differences. The first is the angular 
power spectrum, defined as the ensemble average Ci =< |a/ m | 2 >, where the temperature on 
the sky is expanded in spherical harmonics as (5T/T)(8,<f)) = Ea/ m Y; m (# , <fi) . In Figure 5, 
the quantity /(/ + l)Ci, constant for a perfectly scale invariant spectrum, is plotted against 
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/. Multipole index / is converted into an angular scale by noting that on small scales, the 
sky is nearly flat, and I is effectively the two dimensional wavenumber. Modes with index I 
have a wavelength on the sky of A6 = (2n/l). Figure 5 shows the ensemble averaged power 
spectra for the four theories considered here. Each has been normalised to the COBE 4 
year data, using the results of [EC for monopoles and texture, and |TJ| for strings [JTBJ . The 



curves show the average of 3,4 and 6 runs (producing 8 maps each) for strings, monopoles 



and textures respectively. The fractional statistical error in the power at each I is ±y2/iV 
for iV uncorrelated maps - for the available ensembles this works out at approximately 20 
per cent for textures, 25 per cent for monopoles, 30 per cent for strings. The dashed line 
shows the result for the standard inflationary model. 

The solid line shows the result for the texture model calculation presented in 0, when 
it is used to compute the same contributions to the anisotropy studied here. The CVs from 
the model have been multiplied by a factor 0.6 after COBE normalisation (this correction 
is consistent with the estimated systematic errors quoted in ||) - it appears that model 
overestimates the height of the primary peak by this factor. The position of the primary 
peak is however in excellent agreement with the predictions of refs. |5|], ||. Note that the 
approximation employed here of neglecting the Sachs- Wolfe integral significantly underem- 
phasises the second peak (at I ~ 650), but not the third, which also appears visible in the 
full texture simulations. 

The higher peaks do appear to be present in texture simulations at the locations expected 
from the model of ||, but the statistics so far obtained are inadequate to fully resolve 
them. One should also worry about effects due to the non-scaling initial conditions, and 
other sources of noise which would tend to smear the peaks. The techniques used here are 
not ideal for quantitatively determining the extent of the decoherence of the defect-induced 
perturbations (nor indeed the C/'s themselves), for which a Fourier space approach would be 
preferable. Nevertheless, secondary peaks are apparently visible in the intrinsic temperature 
term alone (dashed lines). There are also hints of secondary peaks in the cosmic string 
theory (the spacing between the spikes, Al ~ 280 is what one expects for coherent acoustic 



oscillations ) ~ 7r, and A/ ~ A(/cto) ~ 50A(/cr rec ) ), but the statistics are 

again insufficient for any firm conclusion. It is notable that the power spectra for individual 
maps frequently show very clear peaks, with the appropriate spacing Al ~ 280, but these 
are displaced from map to map in such a way that the ensemble average is smoother. This 
seems to be a clear signal of decoherence J7J. 

The calculations reported here, as mentioned above, significantly underestimate the total 
power at lower values of /. The scale invariant plateau produced by the Sachs- Wolfe integral 
is at 1(1 + l)Ci ~ 7 x lCT 10 . For 300 < / < 2000 or so for textures and monopoles, and 
400 < / < 2000 for strings, the CVs shown include the dominant effects, but the Sachs- 
Wolfe integral may add ~ 20 — 40 per cent to the power. Above I ~ 2000, at least in the 
case of strings, the structure of the defect-induced Sachs- Wolfe signatures produced after 
last scattering start to contribute, leading to a tail of high I power which is unaffected 
by Silk damping [[1{J. It is interesting to compare results with Bouchet et. al. p0fl , who 



calculated only the Sache- Wolfe integral for strings, in a Minkowski-space approximation. 
They concluded that the contribution was roughly scale invariant, and the rms anisotropy 
from the interval [Z, 2Z] was ~ 6G/i, with G Newton's constant and \i the string tension. 
Here I find the rms from the the intrinsic and Doppler terms to be larger, ~ 13G/i. Thus 
the contributions studied here are likely to mask the linear signatures of the strings, at 
least on angular scales above 10 arc minutes. These conclusions are consistent with model 



calculations of the string anisotropy power spectrum [17] . 

The number of maxima and minima of height ua, where a is the standard deviation, 
provides a simple quantitative test of the nonGaussianity of the maps. For a Gaussian 
theory, the number of peaks and minima per unit area, dN/ dv is straightforwardly calculable 
in terms of moments of the power spectrum |[21|| . In this case the curves for maxima and 
minima are simply related by the replacement v — > —v. Figure 6 shows the number of 
maxima and minima for the theories studied here. While the curves for minima take a form 
very similar to what one expects for a Gaussian theory, the curves for maxima show a clear 
excess of high peaks, especially marked in the monopole and texture theories. Note that 
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this skewness was not present in the maps made by Coulson et. al. |14| of the Sachs- Wolfe 
contribution, it is entirely an effect of the 'intrinsic temperature' and 'Doppler' terms studied 
here. 

One can read off from Figure 6 how much sky coverage is needed in order to see one of 
these extreme events. For example, there is approximately one peak in excess of 5cr in the 
texture theory in each 300 square degrees, and in the monopole theory in each 100 square 
degrees. By contrast, the inflationary theory predicts on average less than one such peak on 
an entire sky. The apparently skew shape of the texture and monopole maxima curves near 
the peak may also be a useful discriminator. 

In conclusion, the prospects for distinguishing between the inflationary and defect theo- 
ries appear excellent. The power spectra of fluctuations are very different, with the defect 
theories having the primary peak in the angular power spectrum shifted to smaller angu- 
lar scales. And the nonGaussianity of the defect maps, with an excess of high peaks on 
scales of order 20 arc minutes, should be readily detectable by the VSA experiment or the 
COBRAS-SAMBA satellite. 
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FIG. 1. A ten degree square patch of the microwave sky according to the cosmic string theory 
for structure formation. The colour scale shows the temperature contrast 5T/T in units of one 
standard deviation. The numbers along the sides of the Figure are grid units. The standard 
deviation for the string maps, adopting the COBE normalisation of ref. [14], but with a sky rms 
[18] on 10° of 35 /j,K, is a = 2.5 x 1(T 5 . 
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Intrinsic Anisotropy from Monopoles 




FIG. 2. As in Figure 1, but for global monopoles. The standard deviation for the monopole 
maps, adopting the COBE normalisation of ref. [11], but with a sky rms [18] on 10° of 35 fiK, is 
a = 2.2 x 1(T 5 . 



13 




14 



Intrinsic An isotropy fro nn Standard CDfvl 





100 1000 1 

FIG. 5. The angular power spectrum of intrinsic CMB anisotropics in the cosmic defect 
theories. The anisotropy is expanded in spherical harmonics, (5T/T) = T,ai m Yi m (9 , <fi) , and Ci is 
defined as < |a/ m | 2 >. The solid red, green and blue lines, linking crosses, triangles and boxes 
respectively, show the numerically determined power spectra for cosmic strings, global monopoles 
and textures, when the theories are normalised to COBE according to the results of refs. [11] and 
[14]. The dashed lines show the contribution from the intrinsic temperature term (5T/T)i alone 
(see text). The solid black curve shows the results for the texture model of reference [5], scaled 
down by a factor of 0.6 after normalising to COBE at low /. The dashed line shows the full Ci 
spectrum for the 'standard' inflationary theory. 
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FIG. 6. The differential number of maxima and minima of the temperature maps in a given 
range of (5T/T) = z/cr, measured in units of the standard deviation a. The solid red, green and 
blue lines, linking crosses, triangles and boxes respectively, show the results for maxima for cosmic 
strings, monopoles and textures. The dashed lines show the number of minima plotted against 
the negative of the height in the same theories. In a Gaussian theory, the curves for maxima 
and minima are identical, and are computable from the power spectra G\ [21]. The curves shown 
for minima resemble those for a Gaussian random field, whilst the curves for maxima are clearly 
nonGaussian, most markedly so for monopoles and textures. 
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